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ABSTRACT 

We present new multi-wavelength observations of the first submillimetre-selected galaxy 
(SMG), SMMJ02399-0136 at z = 2.8. These observations include mapping of the CO 
J=l— i>0 emission using elements of the Expanded Very Large Array, as well as high- 
resolution 1 .4-GHz imaging and optical/infrared (IR) data from the Very Large Array, Hubble 
Space Telescope, Spitzer and Keck-I. Together these new data provide fundamental insights 
into the mass and distribution of stellar, gas and dust distribution within this archetypal SMG. 
The CO J=l— >0 emission, with its minimal excitation and density requirements, traces the 
bulk of the metal-rich molecular gas, and reveals a molecular gas mass of '^lO^^ M0, ex- 
tending over approximately 5 arcsec (^25 kpc in the source plane), although there is tentative 
evidence that it may be significantly larger. Our data suggest that three or more distinct struc- 
tures are encompassed by this molecular gas reservoir, including the broad-absorption-line 
(BAL) quasar from which the redshift of the SMG was initially determined. In particular, the 
new rest-frame near-IR observations identify a massive, obscured, starburst which is coinci- 
dent with a previously known Lya cloud. This starburst dominates the far-IR emission from 
the system and requires a re-assessment of previous claims that the gas reservoir resides in 
a massive, extended disk around the BAL QSO. Instead it appears that SMM J02399-0136 
comprises a merger between a far-IR-luminous, but highly obscured starburst, the BAL QSO 
host and a faint third component. Our findings suggest that this archetypal SMG and its im- 
mediate environment mark a vast and complex galactic nursery and that detailed studies of 
other SMGs are likely to uncover a similarly rich diversity of properties. 
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SMM J02399— 0136 — galaxies: starburst — infrared: galaxies — radio continuum: galaxies 



1 INTRODUCTION 

Submillimetre- (submm-)selected galaxies (SMGs) were discov- 
ered using the Submm Common-User Bolometer Array (SCUBA - 
iHolland et &\\ \\999). expl oiting gravitational amplication by mas- 
sive, foreground clusters jSmail et"al]|l997l l2002t) . They are re- 
sponsible for a significant fraction of the cosmic submm and far- 
IR background jpixsen et al.lll998l) and hence constitute an impor- 
tant class of distant galaxies. The Very Large Array (VLA) facili- 
tated ide ntification of the first SMGs, via high -resolution 1.4-GHz 
imagmg dSmail et alj[200(]| : llvison et aLlbOOd) enabling the acqui- 
sition of spectroscopic redshifts, e.g. z = 2.8 for the fi rst SMG 
to be discovered, SMM J02339-0136 jlvison et al.l 19981 hereafter 



SMM J02399 and 198). The rest-frame UV spectral properties of 
SMM J02399 indicated that th is sy stem contained a low-ionisation, 
BAL quasar, denoted LI by Il98l with a low-surface-brightness 
companion (possibly due to scattered AGN light), L2, associated 
with a larg e Lyg-emitting cloud, extending over at least 20 kpc 
(198; Iviso n et alJl999l:lyillar-Martfn et al.ll999l : IVernet & Cimattil 
1200 li) . iVemet & Cimattil (2001) report an absorption complex at 



— l,000kms~ and, for high-ionisation lines, at — 6,700kms~ , 
both measured relative to LI (at z = 2.795, which is ~ 270 km s~^ 
blueward of L2). Mid-IR and X-ray data have suggested that the 
immense luminosity of this system, Lboi ~ 1-2 x 10^'' L©, 
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arises roughly equally from the AGN and starburst c omponents 
teautz et al.ll20od iLutz et alJlOO^ : IValiante et all2007l) . 

SMMJ02399 w as also the first SM G whose molecular gas 
was detected via CO jpraver et alj|l99^ . This exploited the large 
positive A'-correction afforded by even mildly-excited high- J CO 
lines and the redshifting of t heir frequencies to transmissive parts 
of the atm osphere ( see also Brown & Vanden Boud 1 199 ll , 1 19921 : 
ISolomon et al. 1992: ISolornon & Vanden Bout 200^). Since then, 
great progress has been made in discerning the H2 gas mass dis- 
tribution in distant SMGs, helping to characte rise their evolution- 
ary state jpraver et alJl999l : llvison et al.l200ll) ; recently, more sen- 
sitive and higher resolution work has employed the Institut de 
Radioastronomie Mill imetr ique's (IRAM's) Plateau de B ure In- 
terferometer (PdBI: fPownes & Solomo n 120031: iNeri et al 



iGreveetalJboOSi : iRneib et alj i2005i : .Tacconi et alj 
On the basis of CO J=3— >^2 observations, iGenzel et al 
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hereafter G03, argued for the presence of a massive and rapidly 
rotating disk in SMMJ02399, associated with Ll, although they 
could not rule out an alternative configuration with two galaxies 
orbiting one another. They reported a ~ 3-arcsec structure in their 
sensitive rest-frame 335-^m continuum imaging, associating this 
withLl. 

Only a handful of SMG stu dies dOreve et al.l l2003l : 
iHainline et a"i]|2006l : ICarilli et aljboid) have made use of the CO 
J=l— >-0 line - the transition most capable of tracing the bulk of 
the H2 gas, with its minimal excitation requirements (ncrit(H2) ~ 
410 cm~^ and Eio/ka ~ 5.5 k). Instead, high- J (and thus high- 
excitation) CO transitions have usually been used. CO J=3— )-2 
and higher- J rotational levels probe the star-forming molecular gas 
(e.g. IZhu et aljlioosi : lYao et al. ''2003) and thus give a potentially 
biased view of the total gas mass and gas distribution. 

It must be noted that even when it comes to the proper- 
ties of the star-forming gas, comparisons between SMGs and lo- 
cal ULIRGs are difficult since the high-J CO lines routinely ob- 
tained for SMGs have only recently become av ailable for the global 
molecular gas reservoirs of ULIRGs (e.g. IPapadopoulos et akl 
l2007h . Observations of the CO J=l— >^0 line, on the other hand, 
allow the most straightforward estimate of metal-rich H2 gas mass, 
its distribution, and the enclosed galactic dynamical mass (assum- 
ing dynamically relaxed systems). Comparisons with the available 
CO J=3— )-2 also allow an assessment of the global molecular gas 
excitation conditions, which are in turn a powerful indicator of the 
star-forming and non-star-forming molecular gas fractions in galax- 
ies. The latter has been found to vary c onsiderably, even among 
vigorous starbursts i n the local Universe 1 WeiC et alll2005h and at 
higher redshifts (e.g. IPapadopoulos & Ivisonll2002l) . 

Here, we report sensitive new interferometric observations 
of SMMJ02399, exploiting the arrival of Ka (26.5-40 GHz) Ex- 
panded VLA (EVLA) receivers at the VLA to trace the bulk 
of its molecular gas through its CO J=l— >-0 line emission red- 
shifted into a hitherto inaccessible frequency range (over half of 
radio-identified SMGs lie at z = 1.8 8-3.35 corresponding to CO 
J=l— >-0 in the Ka frequency band, IChapmanetaLll200^ . This 
well-studied and bright target provides an obvious starting point 
for future CO J=l— >-0 surveys of the SMG population with EVLA. 
We also present sensitive, new 1.4-GHz continuum imaging, spec- 
troscopic imaging in the rest-frame UV from the Keck-I telescope, 
and deep multi-filter imaging with Spitzer and the Advanced Cam- 
era for Surveys (ACS) aboard the Hubble Space Telescope. 

Throughout this paper we assume a flat cosmological model 
with D.A = 0.7 and Ho = 70kms"^ Mpc~\ 



2 OBSERVATIONS 

2.1 CO J=l— ^0 observations using the EVLA 

Around 70 hr of data were obtained with the VLA in its C config- 
uration during 2009 July-September, around half via dynamically- 
scheduled 2-hr blocks, the rest via fixed 6-hr blocks. The weather 
conditions were generally excellent. We discarded 20 hr of data, 
taken when a pointing bug was affecting the array, and around 6 hr 
of data, taken when the weather was unsuitable. 

A compact calibrator, 0239-025 (1.0 Jy at 30.3 GHz), lies 
~ 1 deg from SMM J02399, and we used this to track the phase 
and amplitude every 3-5 min, as well as to update the VLA's point- 
ing model every hour (at 8.4 GHz). 0137+331 (3C48) was used to 
calibrate the flux density scale. 

Typically, around ten receivers were functional at any one 
time, each mounted on an upgraded EVLA antenna. We used con- 
tinuum mode, giving dual-polarisation channels, each with a band- 
width of 50 MHz (or ~ 500 km s"^ for CO J=1^0 at ^ ~ 2.805). 
The observed frequency of the CO J=l— >-0 transition lies below 
32 GHz for 2 ~ 2.8, where only two intermediate frequencies 
(IPs), B and D, can be tuned, so we placed IPs A and C at 32.3 GHz 
to constrain the off-line flux density, even though we expected no 
significant continuu m em ission from either the synchrotron, the 
free-free or the dust ( Il98l) . We alternated the tuning of IPs B and 
D be tween the blue and red peaks seen bvi Praver et al. Ul99^ and 
IgoI although some data were taken at non-optimal frequencies. 

The data were processed using standard AXVS procedures 
developed for high-frequency observing. Our calibration strategy 
was able to track changes in phase and amplitude accurately across 
the full 0-300 kA of ui) coverage. Imaging was accomplished using 
IMAGR within AIVS. 



2.2 CO J=3^2 and 1.3-nmi continuum observations 



lG03f s 1.3-mm continuum imaging of SMM J02399 from IRAM's 
PdBI, and their CO J=3— >-2 data cube, were kindly made available 
to us by L. Tacconi and A. Baker (private communication). We use 
these data here to facilitate the most complete possible picture of 
the molecular gas in SMM J02399 and its physical conditions. The 
1.3-mm imaging yields a beam of 1.8-arcsec x 1.4-arcsec (very 
close to that of our 1.4-GHz radio map, see ^2.3\ , while the CO 
J=3— >-2 map has significantly lower resolution, 5.2-arcsec x 2.4- 
arcsec. A full description of the observations and their reduction is 
given in lG03l 



2.3 1.4-GHz imaging with the Very Large Array 

Data at 1.4 GHz were obtained for the Abell370 cluster field dur- 
ing August and September in both 1994 [B configuration] and 
1999 [A], recording 7 x 3.125-MHz channels in each of two dual- 
polarisation IP pairs. 0259+077 [B] and 0323+055 [A] were used 
to track amplitude and phase and 0137+331 was used for flux 
and bandpass calibration. Data processing and imaging followed 
the description given bv lOwen & Morrison (2008). The result- 
ing continuum map has a noise level of ~ 5.7 /iJy beam" ^ near 
SMMJ02399, with a 1.8-arcsec x 1.6-arcsec synthesised beam 
(position angle, PA, 166 deg). 
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Figure 1. True-colour representations of the HST/ACS SM4 ERO and 
NICMOS observations of SMMJ02399 (comprising F475W, F625W and 
F160W exposures) as seen before (inset) and after (main panel) coirecting 
for the gravitational lensing model described in i]2.7l In the rest-frame UV 
imaging, we identify the BAL QSO, LI, as well as a marginally extended 
and somewhat bluer companion, LIN, and the inuch bluer L2. In contrast, 
the NICMOS image identifies an extended source, L2SW, which is very 
red in the rest-frame UV-optical. This source appears to coincide with the 
far-IR-bright and gas-rich emission region within this system. These four 
primary einission regions are marked. 



2.4 Hubble Space Telescope and Spitzer observations 

Following the successful Servicing Mission #4 in 2009 May, ob- 
servations were obtained using the Advanced Camera for Surveys 
(ACS) during 2009 July (PID 1 1507), for a total of 6,780, 2,040 and 
3,840 s in the F475W (g), F625W (r) and F814W (i) bands, respec- 
tively. The data were reduced using the MULTIDRIZZLE package 
and the images are availablj^- 

In addition, we have retrieved Near Infrared Camera and 
Multi-Object Spectrometer (NICMOS) imaging of SMMJ02399 
(PID 11143). These were reduced and combined in the standard 
manner from the pipeline-reduced data in the Space Telescope Eu- 
ropean Coordinating Facility archive. They comprise 2,560 s of in- 
tegration with the NIC2 camera in the Fl lOW ( J) and F160W (H) 

filters. 

Finally, Abell370 was observed wi th Spitzer (IWemer et al.l 

l2004h in all the IR Array Camera (IRAC - iFazio et alj2004h bands 
(3.6-8.0 /xm) in 2006 February (PID: 137). The total integration 
time was 9 ks in each of the four bands: 3.6, 4.5, 5.8 and 8 ^tm. 



2.5 Keck spectroscopic observations 

Keck-I long -slit spectroscop y was obtained during 2001 October 
using LRIS ( IOkeetal.1 19951) to step through seven positions across 
the object, integrating for 2ks in each, with a 1-arcsec slit and the 
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low-resolution 300 lines mm grism. A full description of the ob - 
servations and data reduction are presented in lSantos et al.l ( |2004|) . 

2.6 Astrometric alignment 

The bright point source, LI, visible in many wavebands, has al- 
lowed us to remove the positional ambiguity that would otherwise 
hamper a detailed comparison of the kind undertaken here. 

We adopt the VLA 1.4-GHz map as our fundamental co- 
ordinate frame and shift the other frames to this one to align 
the emission from LI. We employed local phase calibrators (1.0, 
10.6 and 13.0deg from SMMJ02399), with positions known to 
--^1 milliarsec, and residual phase errors after calibration were 
small, so we estimate the absolute alignment of the 1.4-GHz and 
CO J=l— >0 coordinate frames should be better than 0. larcsec, 
leaving aside uncertainties relating to signal-to-noise and beam size 
(discussed in ^3.4\ . The astrometry of the HST images was firstly 
calibrated against a wider-field optical image of the cluster which 
had in turn been aligned onto the FK5 coordinate system defined 
by United States Naval Observatory optical catalogue to an r.m.s. 
precision of ~0.2 arcsec. To align the IRAM data (1.3-mm contin- 
uum and CO J=3-->2) with the 1.4-GHz continuum emission, we 
shifted the former to the east by 0.5 arcsec in R.A., which is consis- 
tent with the claimed astrometric precision of the PdBI images. To 
align the optical/near-IR data with the 1.4-GHz continuum emis- 
sion, we shifted the former in Dec. by -1-0.15 arcsec. 

2.7 The lensing model 

Having described the large collection of new and archival data at 
our disposal, we describe now how we deal with the gravitational 
lensing experienced by SMMJ02399. We correct for the lensing 
distortion and magnification using a mass model of the Abell 370 
cluster. This is discussed in detail by iRichard et all ( [2009,) and is 
based on ten multiply-imaged systems detected in the HST/ACS 
data ( i]2.4| (. We derive a total magnification factor of 2.38 ± 0.08 
at the location of SMM J02399 in a direction roughly aligned with 
L1-L2 and we use the model to reconstuct the m orphology in the 
source plane using LENSTOOL i Kneib et al.l 19961) . The effect of the 
correction can be seen in Fig.[T] 



3 ANALYSIS AND RESULTS 
3.1 Taxonomy 

We begin with a brief description of the SMM J02399 system, as 
seen in the rest-frame UV-optical true-colour imaging from HST. 
In what follows we will refer to four main components seen in 
Fig-E LI - the BAL quasar to the west, the brightest UV emitter in 
the region; LIN - the bright, compact source, ~8 kpc north of LI; 
L2~ the low-surface-brightness feature, ~ 10 kpc east of LI; L2SW, 
the red feature seen in our NICMOS imaging extending ~6 kpc to 
the south-east of LI. 



3.2 Radio continuum morphology 

In Fig. |2] (central panel) we see the 1.4-GHz synchrotron emission 
overlaid as green contours on the HST true-colour image, together 
with 50-kA-tapered CO contours. A knot of synchrotron emission 
is coincident with LI but the majority of the radio emission extends 
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Figure 2. Left: A true-colour 3.6, 4.5 and 5.8-/.tm image of SMM J02399, where we have removed the contribution from LI by scaling and subtracting a 
point source from each of the bands, revealing emission coincident w ith L2SW. We ove rlay on this the 1.4-GHz continuum contours (green) and the positions 
of the 850-/xm and X-ray centroids (red circle - 198; blue square - iBautz et alj|200dl as revised by |GO30. Middle: CO J=l— >0 (white contours) from an 
image made with a Gaussian taper of 50kA, plus 1.4-GHz continuum (green), overlaid on the HST imaging. CO contours: —3, —2, 2, 3, 4, 5, 6 X c, where 
<T = 50 ^lJy beam~^ ; 1.4-GHz: —4, —2, 2, 4, 6, 8, 10 X a, where a = 5.7 /.tjy bearn"^. R ight: Contours (white) of the CO J=3—^2 velocity-integrated 
emission superimposed on a greyscale image of the 1.3-mm continuum emission jG03l) . lG03l suggested that LI is the primary repository of SMM J02399's 
dust and gas. It is clear that the dusty gas reservoir is not centred on LI, but is coincident with L2SW - the very red source highlighted by the NICMOS 
imaging. 
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Figure 3. Keck/LRIS spectroscopy of the four regions identified in Fig. \T\ 
All four regions display strong Lya emission and - with the exception of 
LI - weak or undetected rest-frame UV continuum emission. The spec- 
trum of LI displays the strong, self-absorbed emission lines characteristic 
of BAL QSOs. The three other regions display Lyo and weak N V emis- 
sion, suggesting local ionisation sources may reside in these components. 
The wavelengths of potential emission features are marked. 



to the east, close to the extended near-IR source, L2SW. It is coinci- 
dent to within the uncertainties with the original SCUBA 850-/im 
continuum centroid, with the PdBI 1.3 -nim cont inuum emission, 
and with the CO J=3— >2 line emission as shown in the 

right-hand panel of Fig.|2] 

There are two possible explanations for the radio morphology 
of this system. First, the radio morphology could be due to a bipolar 
jet from LI - a known X-ray-bright BAL quasar - which encoun- 
ters only one working surface: the cloud of molecular gas traced by 



the CO and dust emission, to the east of LI. In this scenario, L2 
and L2SW may be shock-excited regions - L2 lying at the point 
where the jet terminates - wi th the entire system po wered by the 
kinetic energy of the jet (e.g. IPunslv & Tinga\]|2005l) . Taking this 
int erpretation , SMM J02399 is reminiscent of the system described 
by Jo zsa et alJ ( 120091 ) . where extended emission originating at the 
centre of IC 2497 at z = 0.05, probably a radio jet, points in the di- 
rection of Hanny 's Voorwerp - a region in a neighbouring H I cloud, 
first seen because of its shock-excited, optical emission lines. 

Alternatively, and more plausibly, the extended emission may 
be associated with a heavily dust-obscured starburst between LI 
and L2. As we shall see in the next section, the host of this star- 
burst, L2SW, is identified through IRAC and NICMOS rest-frame 
optical/near-IR emission, as well as wispy rest-frame UV emis- 
sion structures representing leakage from patches of low obscura- 
tion (similar to th e scenario proposed for SMM J14011+0252 by 
llvison et"al]|2001l) . The spatial coincidence of this source with the 
CO reservoir, as well as the 850-/xm and 1.3-mm continuum emis- 
sion, strongly suggests it is the seat of the hyperluminous starburst 
within this system. 

As an aside, the lack of emission at 32.3 GHz rules out the 
possibility that the AGN within LI has a flat-spectrum core. 



3.3 Rest-frame UV/optical/near-IR light 

The most remarkable feature of the 2-D spectroscopy described in 
i]2.5| is that Lya is seen across the entire 1,600-kpc^ region shown 
inFig.[T] The scale of the Lya nebular was hinted at by early data 
( Il98l) but we can see now that SMM J02399 is a classic 'Lya blob' 
(e.g. iMatsuda et alj|2004l) . There is a velocity gradient in Lya as 
one moves from LI to L2. Secondary peaks in Lya intensity and 
velocity dispersion are coincident with L2 - the local velocity gra- 
dient is around ~ 500 km s^ ^ . UV continuum is co-spatial with the 
line emission over much of the region, but none is seen near L2. 
Fig. [3] shows 1-D spectra of the four regions identified in 
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Figure 4. Amplitude of the real visibility components as a function of base- 
line length. The amplitude of the imaginary components, whose l-tr range 
is indicated in grey, is consistent with zero. A Gaussian consistent with a 
source size of 10 kA (20arcsec FWHM)is shown to guide the eye. 



Fig. [T] Strong Lya and N V emission domina te th e spectrum of 
LI, along with weak Hell (A164.5nm) (e.g. "198) at 625.17 ± 
0.20 nm. L2 and L2SW are practically indistinguishable, spectro- 
scopically. Regarding the Lya emission, and the proximity of an 
AGN/starburst, SMM J02399 is very similar to the Lya nebula pre- 
sented by Scarlata et al. ( 2009). LIN, L2 and L2SW all show weak 
N V emission suggesting that they contain ionising sources. 

LI dominates the SMMJ02399 system in IRAC wavebands. 
We show in the left-hand panel of Fig.|2]the true-colour IRAC im- 
age after subtraction of a point source at the position of LL It is 
immediately apparent that L2SW emits strongly at rest-frame 0.95- 
1.5 /im. 

Using the optical photometry from'l98l and the flux densities 
of L2 determined from the NICMOS and IRAC imaging ( ^2A[ 
see also Figs [T] and |2l, we find that the spectral energy distribu- 
tion (SED) is well fit with a reddened and very young stellar pop- 
ulation (Ay ~ 1.5; ~10Myr old). From this we derive a best- 
fit, rest-frame, lensing-corrected ff-band luminosity for the galaxy, 
Mh ~ —25.4, equivalent to a stellar mass of ~IO"M0 usme 
a light-to-mass ratio appropriate for a young b urst (Lh/M ~ 3 ), 
with an uncertainty of at least a factor of 3x teorvs et al]|2005h . 
This confirms the visual impression (Fig. [2} that L2 corresponds 
to a massive stellar system. We note that the rest-frame 1.5-2 /im 
SED of L2 is flat, suggesting that an obscured AGN is not dominat- 
ing its near-IR luminosity (cf. Hainline et al. 2010, in preparation). 
We stress that the combined SED of LI and L2 is dominated by 
the AGN emission from LI. Only the combination of gravitational 
lensing and the high-resolution imaging used here allows us to de- 
couple the AGN activity from that of the luminous, obscured star- 
burst within this system. Similarly detailed studies of other SMGs 
may also un cover a rich variety in their properties (e.g. llvison et al.l 
I200dl2008l) . 

3.4 Molecular gas 

The other major baryonic component in the system is the cold gas. 
The CO ,7=1— >-0 images mentioned in i]2.1l were made using a va- 
riety of Gaussian uv tapers (125, 75, 50 kA) resulting in synthe- 
sised beam sizes of ~1.7, 2.4 and 3.2arcsec and noise levels of 
26-50 fj.}y beam~^ . The noise level rises as the longer baselines are 



weighted down, but the signal due to CO line emission increases 
markedly. This is best illustrated by the amplitudes of the real and 
imaginary visibility components (Fig.|4} which suggest that the CO 
emission is resolved on baselines longer than ~10kA, with a total 
flux density of approximately 500 /iJy. 

The position of the CO J=l-^0 emission (02''39™51.96 ± 
0."03, -01°36'59.8 ± 0.5" J2000; see Fig.© is consistent with 
the emission seen by SCUBA and PdBI. We cannot rule out a con- 
tribution to the CO J=I— >-0 flux from the quasar, but the majority 
of the emission appears to be coincident with L2SW. 

The peak and total flux densities in the most severely tapered 
CO J=I-5>0 map are 325±50 ^ijy beam"^ and 540± 150 ^iJy. The 
latter value is sensitive to the manner in which the measurement is 
made; however, Gaussian fits and a simple integral over the emis- 
sion region in the dirty image with TVSTAT give results consistent 
with Fig.m We adopt 540 ^Jy as the total CO J=I-^0 flux density 
hereafter. Gaussian fits to the dirty image provide a very rough indi- 
cation of the size of the emission region: 2.0 x 5.0arcsec, with the 
major axis at PA, 150 ± 20deg, implying a lensing-corrected size 
of ~25 kpc. This is only barely consistent with the uv visibilities 
(which suggest a size closer to lOarcsec). Without data from the 
VLA in its D configuration, we have only loose constraints on the 
size of the CO J=I— >0 emitting region and may be underestimating 
the total flux density. 

To correct for t he finite throughpu t of our IFs, we adopt and 
adapt the strategy o f lGreveetalJ ( l2003l) . correcting upwards by 45 
per cent to account for the breadth of the line seen by G03 relative 
to our total bandwidth (~90 MHz, taking into account the bandpass 
rollover), which takes the total flux in CO J=I— >-0 to Si-o dV = 
0.70±0.18Jykms-\ 

As expected, no signal is apparent in the 'off-line' IF pair, 
down to a 2-a limit of 5'32.3GHz < 60 /iJy. This situation does not 
change when tapering the uv data, although the noise doubles when 
a 50-kA taper is applied. 

The velocity/area-integrated CO brightness temperature in the 
source reference frame, Li_o = j^y Tb dAdV, can be ex- 
pressed as 



3.25 X 10" f f ^ -1 2 

Li-o = —— — r— / ^i-odVKkms pc , (I) 

(1 + 2) \vcoJ J 

where Dl is the luminosity distance (Mpc), uco is the rest-frame 
frequency (GHz) of the CO J=I— >0 transi tion and the velocity - 
integrated flux density is in Jykms^^ fe.g. ISolomon et al]|l997h . 
The corresponding H2 mass, A/(H2), is given by XqoLi-o, 
where Xco is the CO-H2 luminosity-mass conversion fac- 
tor. For quie scent environtnents, Xco ^ 5 Ma (K kms~ 
pc^)"^ (e.g. ISoIomon et all 1 19871 ; ISolomon & BarrettI Il99ll) . 
In extreme starburst environments, like those prevailing in 
SMM J02399, Xco ~ 1 Mp (K kms~^ pc^)"^ may be more suit- 
able ( Downes & Solomonll 1998b . which corrects semi-empirically 
for the fact that the molecular gas reservoir in the dissipative 
mergers that typically characterise such starbursts is a continuous 
medium with a significant stellar mass component rather than an 
ensemble of virialised clouds (so that CO line widths no longer re- 
flect only the mass of gas). In this case, M(H2) ~ (1.0 ± 0.3) x 
10^^ M0 (after correcting for lensing amplification). 

An absolute lower limit for the H2 gas mass can be obtained 
by assuming that the ^^CO J=l— >-0 emission is optically thin. For 
local thermodynamic equilibrium. 
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M(H2) 

Ll-Q 



0.08 

f [CO/H2 



Z \ J{T^) 



\ 10- 



M 







K km s ^ pc^ 



(2) 



where To = Ei/ks ~ 5.5 k, J{T) = To{e'°'' - Tbg = 

(1 + 2:) TcMB ~ 10 K (the temperature of the cosmic microwave 
background at z = 2.8), gi — 3 - the degeneracy of level n — 1, 
Z ~ 2 (Tk/To) ~ the partition function, and [CO/H2I ~ 10"* for 
a Solar metallicity environment (Bryant & Scoville 1996). For typ- 
ical star-forming gas, where Tk ~ 40 — 60 K, this yields (Xcg") ~ 
0.65 Mq (Kkms-i pc^)"\ and thus M(H2)min ~ 6 x 10^° Mq 
is the minimum plausible molecular gas mass in this system. Of 
course, CO traces metal-rich gas and in metal-poor environments 
(and/or intense far-UV radiation fields) more molecular gas may 
be present where CO has been dissociated. In this sense, and be- 
cause we cannot rule out that we have missed a more extensive CO 
J=l— >-0 component due to the lack of D-configuration EVLA data, 
the latter value is a firm lower limit. 

The availability of both CO J=l— j>0 and J=3— j>2 allows a 
rudimentary probe of the average molecular gas excitation state. 
This can provide powerful clues about whether a galaxy-wide star- 
burst occurs over most of the available molecular gas reservoir 
(pushing the global r3_2 ratio up), or whether a smaller gas mass 
fraction is involved in fueling the starburst. High global high- 
J/low-J CO line ratios, indicative of galaxy-wide starbursts in- 
volving most of the ambient gas reservoir, have been observed 
in va rious hi^h-redshift systems (e.g. £apadopoulos et al. 2 00^ : 
IWeTB et all2007h . but lo w ratios in vigorously starbursting systems 
have also been detect ed /Papa dopoulos & Ivisonll2002l : lDaddi et al.l 
l2009l : lDannerbauere t al. 2009). 

For the reported CO J=3— >-2 velocity-int egrated line flux of 
I3-2 = / S3-2dV = 3.0 ± O.SJykms"^ ( jpraver et alj|l99g . 
lG03h . the globally-average CO (3-^2)/(l->0) brightness temper- 
ature ratio would be r3_2/i_o = {ux-o/u3-2)^{h-2/h-o) = 
0.48 ± 0.13. This line ratio is well within the domain expected for 
globally star-forming molecular gas, although it lies towards the 
low end of the expected range. Quiescent Milky Way clouds have 
r3_2/i-o ~ 0.2 — 0.3 whilst th e typical average value for star- 
burst nuclei is r3_2/i_o = 0.65 jPevereux et al.|[r994l : I Yao et al.l 
l2003h . A grid of radiative transfer large-velocity gradient (LVG) 
models, constrained by the observed ratio and Tk ^ Tdust (since 
gas is always warmer than the concomitant dust where the gas 
is heated by far-UV ra diati on and/or turbulence), with the latter 
being Tkust > 40k (^198^ yields n(H2) ~300-10^ cm"^, for 
Tit = 40 — 110 k, for various values of the average gas velocity 
gradient dV/ dr. Given that other CO line ratios are not available to 
us, there is an irreducible n(H2) — Tk degeneracy (the higher the 
temperature, the lower the density reproducing a given CO ratio, 
typically). If we further constrain the LVG solutions by allowing 
only those corresponding closely to virialised molecular gas veloc- 
ity fields where 



(dV/dr), 
(dV/dr) 



1.54 



[CO/H2] 
\/a Ago 



{n(H2)) 



1 -1/2 



103 cm- 



(3) 

where K^ir 2> 1 and A'vir < 1 correspond to unbound gas 
motions and unphysical sub-virial motions, respectively, we find 
that the solutions with Tk ~ 45 — 50 K are mostly favoured. We 
note, however, that super-virial solutions are plausible if the CO 



J=l— >-0 emission originates predominantly from a cold, extended 
gas reservoir; indeed, K^ii ^ 1 solutions occupy the majority of 
the n(H2) — Tk parameter space compatible with the observed line 
ratio, typically with n(H2) ~ 10^ cm"^ and Tk ~ 55 - 90 K, con- 
dition s that often characterise the gas in starbursts (e.g. lAalto et al.l 
Il995h . 



4 DISCUSSION AND CONCLUDING REMARKS 

Based on the rest-frame UV and optical imaging, it appears that the 
SMMJ02399 system comprises at least three significant compo- 
nents: LI, LIN and L2SW. Are these all massive baryonic galaxies 
within a very large reservoir of molecular gas, each caught during 
a different phase of their development? Or does LI dominate the 
system, with L2SW, LIN (and L2) representing mere outflows and 
reflections? 

G03 argued that the dust and gas seen towards SMM J02399 
are associated with LI and that the dynamics of the system as a 
whole are best described as a massive disk. Based on the data de- 
scribed here and elsewhere, there is little doubt that LI does rep- 
resent a significant concentration of mass: it is the brightest rest- 
frame UV emitter in the region; it is bright in the rest-frame 0.95- 
2. 1 ^m waveband a nd ho sts an X-ray-bright BAL quasar. However, 
Fig. |2] suggests that lG03l 's main assertion may not correct because 
the mm/submm, radio and CO emission lie predominantly on and 
around L2SW, although some weak emission does appear to be as- 
sociated with LI . 

Moving on to L2 and L2SW, IVemet & Cimattil ( l200lh and 

GO3! argued that L2 is due to scattered light from LI, although 
G03 did concede that the L1/L2 system may comprise two galax- 
ies orbiting one another. As we have argued, the true-colour Spitzer 
imaging presented in Fig.|2]shows that L2SW is a luminous source 
at rest-frame 0.95-1.5 jim. We stress that several of these bands 
are not expected to suffer strong contamination by emission lines, 
hence we are seeing continuum emission, not scattered light. 

We conclude that SMM J02399 comprises at least two galax- 
ies - possibly as many as four - caught at different po ints during 
their e volution. In the evolutionary scheme outlined bv lPage et al.l 
1 I2OO4I) , LIN may be the most evolved of these systems, followed 
by LI, L2 and the gas- and dust-rich system, L2SW, whose rest- 
frame far-IR emission led to the system's discovery. However, this 
sequence could just as easily be related to mass: we have either one 
or two UV systems with modest stell ar mass, the QSO host g alaxy, 
which must be massive according to Magorrian et al.l ( Il998l) . plus 
a massive, obscured starburst. 

The major result from our study is the discovery of a mas- 
sive, cold and apparently extended gas reservoir in SMMJ02399. 
Tentative evidence for luminous CO J=l— »0 emi s sion in SMGs 
had b een publishe d previously ( Greve et al. 2003|; iHainline et al.l 
good) . lGreveetal.1 used the VLA to observe ERO J16452-I-4626.4, 
a submm-bright starburst at z = 1.4, findi ng r2-i/i-o = 0.6± 0.2 
and r5_4/i_o = 0.10 ± 0.05. Similarlv. Mainline et alj observed 
the AGN-dominated z = 3.4 SMG, SMM J13 120-1-4242, using the 
Green Bank Telescope, and found a brightness temperature ratio, 
r4_3/i_o = 0.26 ± 0.06, indicative of a massive, sub-thermally 
excited molecular gas reservoir. Our new observations suggest that 
massive, cold gas components must be common amongst the SMG 
population. 

In retrospect this discovery shouldn't come as a great surprise. 
As is widely appreciated, optical/IR lines can lead to misleading 
views of velocity fields and galaxy types/sizes. They trace ionised 
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gas, whose distribution traces the star formation and AGN activ- 
ity within the galaxy, rather than its underlying structure. Equally, 
however, the excitation requirements of high- J CO lines tie them 
to a particular star-formation area, so they do not give a faithful 
picture of the true H2 mass and dynamical mass distributions in 
galaxies either. The compact ( ^0.5-1 kpc) emitting regions thereby 
revealed l lTacconi et ai]|2008h are unlikely to trace the true dimen- 
sions of galaxies. It can be argued that high- J observations simply 
return the star-formation efficiency (SFE) per star-forming H2 gas 
mass, rather than the more fundamental SFE per total H2 gas mass. 
Indeed, previous high- J CO studies have uncovered low r5„5/3_2 
and r7_6/3 -2 brightness temperature ratios (~0. 15-0.35) in sev- 
eral SMGs jXacconi et al.l 2006) and such a low global high- J/low- 
J CO ratio is exactly how a massive, low-excitation H2 gas phase 
would betray itself at high redshifts. 

A similar situation exists at z ~ where low-excitation H2 
gas is found to ex tend well beyond compact nuclear starbursts (e.g. 
IWeiB et alj|2003) while containing the bulk of the H2 gas mass. 
Line ratios comparable to those seen in the SMGs can be recov- 
ered by degrading similar observ ations of local starbursts and vig- 
orously star-forming spirals (see IPapadopoulos & Seaauistlll998l) 
to the spatial resolution typically available in CO observations of 
the distant Universe (HO kpc). 

The confirmation of the presence of significant amounts of 
cold gas in SMGs has several significant implications for the inter- 
pretation of this population. In particular, the large gas reservoir de- 
tected in CO J=l— >0 highlights the danger of drawing far-reaching 
conclusions from more biased high- J studies, particularly concern- 
ing the gas mass and gas fraction, the dynamics of SMGs, the 
evolutionary links between SMGs and their likely descendants at 
2: ~ 0, their true star-forming potential (i.e. their potential present- 
day mass), the possible durations of star-forming episodes, and the 
state of their velocity fields (ordered and/or disk-like versus the 
more chaotic fields typical of mergers) - all of these issues are im- 
portant, yet the answers remain largely unknown. 

In the fullness of time, we expect EVLA, the Atacama Large 
Millimetre/submm Array and the Herschel Space Observatory to 
fully compare the excitation properties of gas reservoirs in SMGs, 
ULIRGs and QSOs and the relative distributions of star-forming 
and quiescent molecular gas, relating this to the activity in these 
systems, and their evolu tion, thereby constrain ing theoretical mod- 
els of galaxy formation jSwinbank et afluOOSl) . 
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